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Abstract 

This paper focuses on the geometric phase of general mixed states under unitary 
evolution. Here we analyze both non-degenerate as well as degenerate states. Start- 
ing with the non-degenerate case, we show that the usual procedure of subtracting 
the dynamical phase from the total phase to yield the geometric phase for pure 
states, does not hold for mixed states. To this end, we furnish an expression for the 
geometric phase that is gauge invariant. The parallelity conditions are shown to be 
easily derivable from this expression. We also extend our formalism to states that 
exhibit degeneracies. Here with the holonomy taking on a non-abelian character, 
we provide an expression for the geometric phase that is manifestly gauge invariant. 
As in the case of the non-degenerate case, the form also displays the parallelity 
conditions clearly. Finally, we furnish explicit examples of the geometric phases for 
both the non-degenerate as well as degenerate mixed states. 
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1 Introduction 



The notion of geometric phases was first addressed by Pancharatnam [1] in his pi- 
oneering effort to compare the phases of two beams of polarized light. Although, 
his treatment was essentially classical, the notion of geometric phases was later 
shown to have important consequences for quantum systems. Indeed, Berry [2], 
demonstrated that quantum-mechanical systems could also acquire phases that are 
geometric in nature. He showed that, besides the usual dynamical phase, an ad- 
ditional phase that was related to the geometry of the state space was generated 
during an adiabatic evolution. Soon after, this phase was shown to be precisely the 
holonomy in a line bundle, establishing an even stronger link to the geometry of the 
state space. Essentially, Simon [3] furnished an elegant geometrical interpretation 
of Berry's phase in which he considered a line bundle L over the space of parameters 
M of the system. More precisely, it was shown that L, if endowed with a particular 
connection, can reproduce the Berry's phase when closed loops in M are lifted with 
respect to this connection. 

It was subsequently shown that the adiabatic condition could be relaxed by 
considering the holonomy of a U(l) bundle over the projective Hilbert space of 
the underlying system. By appealing to the notion of gauge invariance, Aharanov 
and Anandan [4] showed that one could define a phase that was the same for all 
gauge-equivalent paths. Gauge-equivalent paths, here refer to the infinite number 
of possible evolutions (curves in the Hilbert space of the system) that project to a 
given closed curve in the projective Hilbert space of the system. 

The abelian nature of the holonomy was also shown to be extendible to the 
non-abelian ones. Wilczek and Zee [5] showed that for Hamiltonians that exhibited 
degeneracies, the adiabatic evolution admitted holonomies that were generally non- 
abelian. The non-adiabatic generalization of this was subsequently furnished by 
Anandan [6]. 
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It is interesting to note that the notion of geometric phases has taken on an 
important role in the area of quantum computation [7]. Indeed, it has been demon- 
strated that one can implement quantum logic gates based solely on the concept 
of geometric phase. The geometrical nature of the phase makes such a realization 
fault-tolerant and thus more robust towards noise. The application has also been 
proposed within the context of non-abelian holonomies. Termed Holonomic Quan- 
tum Computation, the program provides an all-geometrical approach to quantum 
information processing [8]. 

While many of the propositions in this area has been centered around pure states, 
the need to address the issue of geometric phase for mixed states is rapidly gaining 
prominence. Systems that are interacting with other systems, the environment for 
instance, are generally described by mixed states. Again, the strong focus in this 
area is fueled by the promise of realizing quantum logic gates under more realistic 
physical conditions. 

Now, among the first to introduce the notion of geometric phase for mixed 
states was Uhlmann [9] . By considering a purification and the notion of parallelity, 
he furnished a definition for the geometric phase for mixed quantal states. Recently 
Sjoqvist et.al. [10] have introduced a new formalism that defines the geometric 
phase in the context of quantum interferometry. Using the standard setup of the 
Mach-Zender Interferometer, they provided a procedure for calculating the geomet- 
ric phase of mixed state undergoing unitary evolutions. They have also furnished 
the necessary and sufficient conditions for parallel transporting mixed states under 
such evolutions. 

Recently, it was pointed out in Ref.[ll] that the two approaches, in general, are 
not equivalent . More recently, Ericcson et.al. [12] have shown that the conditions 
of parallelity used in the two approaches lead to generically distinct phase holonomy 
effects for entangled systems undergoing certain local unitary transformations. 

In this paper, we consider the geometric phase of general mixed states under uni- 
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tary evolution. Here we consider both non-degenerate as well as degenerate states. 
Starting with the non-degenerate case, we show that the removal of the dynamical 
phase from the total phase is not a trivial subtraction as is the case for pure states. 
In particular, we show that such a procedure is generally not gauge invariant. We 
then proceed to furnish an expression that is manifestly gauge invariant. The par- 
allelity conditions are shown to be easily derivable from this expression. We also 
extend our formalism to states that exhibit degeneracies. Parallelity conditions for 
such states are also obtained. Finally, we furnish explicit examples of the geometric 
phases for both the non-degenerate as well as degenerate mixed states. 

2 The Non-adiabatic Geometric Phase 

In this section we begin by briefly reviewing some basic facts about the non-adiabatic 
geometric phase. Here, we will concentrate on the abelian case, starting with pure 
states ([4], [13], [14]). We begin by considering a quantum system S with states 
belonging to the Hilbert space H. The space of physical states corresponds to the 
non-zero vectors in 7i that are rendered equivalent under the relation 

|V>i >~ |V>2 > iff \i>2 >= e ie \4>i > . (1) 

where 9 is any real number. Technically, for an n-state system with dim(H) =n this 
space is simply a complex projective space of dimension (n — 1): 

V W = = CP n - 1 (2) 

In the language of fiber bundles, the construct of interest is a principal bundle 
over V(H) with a fiber that is isomorphic to U{\). Here the bundle space P corre- 
sponds to the space of normalized state vectors H — {0} with a natural projection 

7r:P^V(H) = M realized through tt : \ip \ip >< (3) 
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Now the unitary evolution of the state vector 

\m>^\m>=w)\m> w 

produces a path C : [0, r] — ► P in the bundle space. This is dictated by the 
Hamiltonian via the Schrodinger equation. The corresponding path Co in the space 
of physical states is obtained by projecting the path C through the map it, i.e. 
7r(C) = Co- Moreover, we say that the evolution is cyclic if the path in M. is closed; 
i.e. if |Y>(t) >< -0( r )l = IV'(O) >< I* i s important to note that there are an 

infinite number of paths in P that project to the same path in M. Indeed, if C and 
C are given by \tp(t) > and \tp'(t) >= e l f^\ij)(t) > respectively for any arbitrary 
real function f(t), then it is easy to see that they define the same path Co in M 
under the projection it. For a cyclic evolution, the total phase acquired corresponds 
to an element of the fiber i.e. e l ^ T £ U(l). The total phase 4>t, here, being the 
argument of < i/j(0)\i/)(t) >: 

total phase 4>t[C] = arg{< ip(0)\ip(T) >} 

= arg{< mMt)\m >} (5) 

In general, this depends on the path C taken in the bundle space P, which in turn, is 
defined by U{t). This is reminiscent of the Pancharatnam analysis which compares 
the relative phase between two states. Aharonov and Anandan [4] showed that if a 
quantity (called the dynamical phase), 

Mc] = - [ T <m\H\m >dt 

Jo 

= -i T < ^(0)|W(t) t W(t)|V'(0) > dt (6) 

JO 

is subtracted from the total phase (5) then the resulting term 

(Pg = arg{< V(0)|W(*)|V(0) >} + i T < ^(0)|W(t) t W(t)|V'(0) > dt (7) 

Jo 

turns out to be a functional of Co only. This is the notion of gauge invariance as 
emphasized in Ref [13]. Since the path C in P from a given intial state |^(0) > is 
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solely determined by U{t) we can label the path by the latter. Here the statement 
of gauge invariance for the geometric phase for all paths differing by a U(l) phase 
reads as:: 

Mmy m ] = ct> G [u{t)} (8) 

for any arbitrary real function f(t). Thus equation (7) furnishes a gauge- invariant 
expression for the geometric phase of a pure state. It is worth noting that arbi- 
trariness of f(t) allows one to impose a condition known as the parallel transport 
condition: 

< ip(0)\U(tfU(t)\ip(0) >= (9) 

which renders the dynamical phase zero. Under this condition the geometric phase 
is just the total phase. 

For mixed states, Sjoqvist et. al. [10] showed that the total phase acquired 
by the system under unitary evolution is a natural generalization of the pure state 
case. Indeed, by considering a Mach-Zender interferometer, in which the interfering 
particles carry an additional degree of freedom (for which the state is mixed), they 
showed that under a unitary evolution 

P (0)^P(t)=U(t)p(0)U(t^ (10) 

the interference pattern of the intensity profile takes the form of 

/ oc 1 + |Tr(W(r)p(0))| cos[ X - arg{Tr(W(r)p(0))}]. (11) 

Here \ is a variable relative U(l) phase in one of the interferometer beams. From 
the shift in the interferometric pattern, they surmised that the total phase acquired 
by a mixed state under unitary evolution is 

7T = arg{Tr(W(r)p(0))}. (12) 



6 



In a basis in which the density matrix p(0) is expreseed as, 

N 

p(0) = J2^k\k><k\, dim(H) = N (13) 
k=i 

the total phase is just the weighted sum of the pure state phases: 

7 t = arg{Tr(^(r)p(0))} 

N 

= arg{]>> fc <k\U(r)\k>}. (14) 
k=i 

In Ref. [10], the parallel transport condition for non-degenerate state was imposed 
by requiring that Tr[p(t)U(t + dt)U{t)^\ be real and positive which in turn leads to 

Tr[p(t)U(t)U(tf] = Ti{p{0)U(tfU(t)} = (15) 

or equivalently 

N 

Y,^k< k\U{tfU{t)\k >= 0. (16) 

k=l 

Under this condition the dynamical phase 

ID = ~ [ T dtTr(p{t)H(t)) 
Jo 

= -i dtTr(p(0)U{tfU{t)) (17) 
Jo 

vanishes identically. It is important to note that if condition (16) is not met then 
the dynamical phase term in eqn.(17) cannot be removed from the total phase term, 
(14) by a simple subtraction; since the latter is defined by a weighted sum of phase 
factors while the former is a weighted sum of phases. The authors of Ref. [10] have 
also asserted that conditions (15) or (16) while necessary, is not sufficient. Instead 
they proposed stronger conditions 

< k\U{tfU{t)\k >=0 Vfc = l,2,..JV (18) 

in which all the constituent pure states in the mixture are required to be parallel 
transported independently. 



In the following section we elucidate the nature of these conditions. In particular, 
we furnish these by first constructing an expression for the geometric phase that 
is manifestly gauge invariant. This is first done for the non-degenerate case in 
which we show that the removal of the dynamical phase requires the use of the 
stronger conditions noted above. We also extend our analysis to states that are 
degenerate. Here, we also provide a generalized version of the geometric phase that 
is applicable when the state is multiply degenerate. The corresponding parallel 
transport conditions are also furnished 

3 Gauge-Invariance and Holonomy in Mixed 
States 

To address the issue of gauge invariance in the context of holonomy for mixed states, 
it is instructive to consider the approach of Boya et.al.[15]. Here one considers 
the little group of a particular mixed state density matrix which can be taken as 
the gauge group for the system that is undergoing unitary evolution. It should be 
noted that the space of states accessible through unitary evolutions from a reference 
state does not cover the entire state space. For instance, if we are considering 
a spin-^ system, the state space consists of all points on and inside the Bloch 
sphere. However, if we start with a particular mixed state and assume that it is 
undergoing unitary evolution then the accessible states lie on a spherical shell with 
radius equal to the magnitude of the initial Bloch vector. In determining the relevant 
bundle structure, we appeal to an important theorem for spaces with transitive 
group action. Essentially, it is well known that if a group G acts transitively on a 
space M then for each x G M there is a bijection between G/G x and M where G x 
denotes the little group or the stability group at the point x [16]. For the problem 
at hand, we are interested in the little group associated with the density matrix at 
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time t = 0. The theorem basically assures us that the orbit space under unitary 
evolution is isomorphic to the coset space, G/G x . More importantly, if the little 
group is closed then (G, G x , ir, G/G x ), where -it : G — > G/G x is the canonical 
map which sends each element of the group to the cosets, serve to define a bundle 
structure. Here G is taken as the bundle space, G/G x the base space, G x the fiber 
and 7r the projection. 

To carry the analysis further, it is instructive to consider a density matrix in its 
diagonal basis (13). Let us start with the case in which the iV-state density matrix 
has eigenvalues that are distinct, i.e. non-degenerate. Then the little group in this 
case is simply 

N 

G x = U(l) x U(l) • • • x U(l) = V e i6 "\n >< n\ (19) 

' N ' ^ 

where {0 n }n=i,2,...N are arbitrary parameters or phases. The group G in this context 
is U(N). It is easy to see that under 

N 

U{t) G U(N) U'(t) = U{t) e i6n(t) \n >< n\ (20) 

n=l 

the orbit of the density matrix remains unchanged: 

p(0) ^ p'(t) = U'(t)p(0)U'(rf = U(t)p(0)U(tf = p(t). (21) 

In other words, there are infinite number of orbits in G that correspond to same 
path for p(t). Now, unlike the pure state case (see (7)), the removal of the dynamical 
phase from the total phase does not render the resulting functional gauge invariant. 
Indeed under a gauge transformation (20) the total phase (14) transforms as 

7T ^ 7 ' T = axg{Tr[p(0)Z/(r)]} 

= axg{5> fc < fc|W(r)|fc > e**W} (22) 
k 

while the dynamical phase (17) transforms as 

ID - l' D = ~i f dtTr(p(0)U'{tfU'(t)) 
Jo 
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= -i [ T dt Tr(p(0)U(tfU(t)) + $> fe fc (r). (23) 
Jo k 

It is evident from (22) and (23) that the ^-dependence in total phase cannot be 

removed by simply subtracting the dynamical phase term. It is also easy to see that 

only in the pure-state case, with only one non-zero u, do we have gauge-invar iance. 

To alleviate this problem we propose a functional, of the following form: 

1G [U] = arg{]>> fe < k\U(r)\k > e ~ £ *<*mm)\k> ]} (24 ) 
k 

which is manifestly gauge invariant (as in (7)). Indeed, the 6 term generated in (22) 
under a gauge transformation (20) is exactly cancelled by a term arising from the 
exponential term in (24). It is easy to check that the above expression corresponds to 
the geometric phase in the case when the state is pure. Gauge invariance also assures 
us that the expression depends only on the path in the state space. This expression, 
thus furnishes us with a gauge-invariant definition for the geometric phase of mixed 
states that are non-degenerate. It is worth noting that when the stronger conditions 
(18) are imposed the resulting expression is nothing but the original total phase 
term, which means, they are justifiably parallel transport conditions; since only 
under parallel transport is the geometric phase the same as the total phase. 

The above analysis holds only for mixed states that are non-degenerate. For 
degenerate mixed states, it is well known that the holonomy group can take on a 
non-abelian character [15]. To this end, let us consider a density matrix that has 
the following diagonal form: 

N 

p(0) = J2^k\k><k\ (25) 
k=l 

in which n of the eigenvalues are degenerate i.e. oj\ = uj 2 = ■ ■ ■ = 0J n with the 
rest distinct. Here the Hilbert space H can be regarded as the direct sum of two 
subspaces, 7i n © H m of dimensions n and m = N — n. The little group of p(0) is 
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non-abelian and has the following form: 



G x = U{n) x U(l) x U(l) x • • • x U(l) 

S v ' 

m 

or in the {\k >}k=i,2,...N basis it takes the form 

a(t) 



(26) 



V(t) 



e i/3 n+1 (t) 



(27) 



Here the only non zero terms are the elements of a, an n x n unitary submatrix and 
exponential of the /3's along the diagonal (m x m)-submatrix. It is easy to see that 
under a gauge transformation 



U{t) ->U'(t) = U{t)V{t) 



(28) 



the path taken by the density matrix remains unchanged. It is worth noting that, 
since both U(t) and U'{t) are the time evolution unitary operators for the system, 
we must require that U(0) = W(0) = 1. As a consequence, the element of the little 
group V(t) must also satisfy V(0) = 1. Now, consider the total phase term (14) 
under such a transformation: 



= aig {Tv n [p(0)U'(r)]} 

= arg{]>> fc < k\U(r)V(r)\k > . 

k 



(29) 



Following the non-denegerate case, we note that in constructing a gauge-invariant 
functional an expression of the form 



7 [W] = arg{Tr w (p(OMrW;T])} 
with a functional T%A\ r] that transforms as 



(30) 



(31) 
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under a gauge transformation would be required. To identify this functional, we first 
observe that in the non-degenerate case, the corresponding functional e~ fo d * <fe l w (*) tz ^(*)l fe> 
can be obtained by fixing the 6 terms in a general evolution matrix (20) through 
the conditions (18) for U' . This yields 



9 k (t) = i f dt' < k\U{t') ] U{t')\k > k = 1, 2, N 
Jo 



(32) 



which upon substitution into (20) leads to the required form. In the degenerate case, 
we have n 2 + m arbitrary parameters that need to be fixed; n 2 for the a matrix and 
m parameters corresponding to the /3's. We can do this by imposing 

<n\U*U'\v> = n,v = 1,2, ...n 
<k\U'^U'\k> = k = n+l,n + 2,...N 

where we have defined basis vectors {\fx >}^=i,2,...n an d {\k >}k=n+i,n+2,...N to 
designate subspaces 7i n and TL m respectively. With 



U' = U{t) 



a(t) 



ij3 N {t) 



(34) 



conditions (33) lead to the following equations: 



$ k (t) = i{U{tfti{t)) kk 



(35) 
(36) 



where are the elements of the a matrix. Equations (35) and (36) can be formally 
integrated to yield 

ft 

aw[W;t] = <fi\Pexp- dt'U{t'yU{t')\v > fi, v = 1, 2, ...n (37) 

Jo 

Pk[U;t] = if dt' <k\U{t') ] U{t')\k> k = n+l,n + 2,...N (38) 
Jo 
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where P denotes path ordering. We can now choose 

a[U;r] 



( e 



i(3 n+1 [U;T] 



(39) 



e if3 N \U;r 

It is worth noting that T\IA\ r] is block diagonal in the chosen basis and that it can 
be written as direct sum form: 



T\U-T\=T Un %A-T\®T Um %A-T\ 



(40) 



with ^„[W;r] = a[U;T] and T h JU;t\ = dmg{e i ^+^ T \e i ^+ 2 ^ T \ ....e^^). 
We next demonstrate that expression (30) is gauge invariant. Indeed, under 



U^U = U(t) 



Vi(t) 



( gi9n+l(t) 



(41) 



g*fjv(t) 

where V\(t) is an arbitrary n x n unitary matrix, we note that 

ThJA;t] - diag(e^+ 1 W,e^+ 2 W,..., e ^W)t^ m [^;r] 



(42) 
(43) 



which implies that j[U] — > 7[W] = t[(W)] (see Appendix for details). Thus expres- 
sion (30) furnishes a natural generalization of the geometric phase for a degenerate 
mixed state. It is easy to see that the conditions (33) here serve as the parallelity 
conditions. 

It is interesting to note that the above analysis can be extended to density 
matrices that are multiply degenerate. For instance if p(0) has eigenspaces Hi, H2 
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...7i m with degeneracies ni, n<i ,...n m respectively, then the Hilbert space of the 
system can be written as a direct sum of degenerate subspaces: 

TO 

H N = Hi®H 2 @ ...®H m with N = J2 n k- (44) 

k=i 

The corresponding T- functional for Hn can also be written as the direct sum of 
.F-functionals define on the subspaces: 

T Un [U, t] = F Hl [U, r] © T U2 [U, r], ©, ®T Um [U, r] (45) 

where the components of jF^ fe [W,r] are given by 

Th&Apm =< Mfe|P exp- [ T dtU{tfU{t)\v k > fi k ,v k = l,2,...,dim(W fc )(46) 

Jo 

if dim (Wfc) / 1 and 

otherwise. The geometric phase is then calculated directly from 

7 [W] = a rg{Tr HN (p(0)U(T)f HN [U;T])} (48) 

4 Explicit Examples 

In this section we furnish explicit examples of the geometric phase for both the 
non-degenerate as well as the degenerate mixed states. We begin with the non- 
degenerate case. To this end, we consider the simplest case of a spin- \ system with 
a density matrix given by 

p= l -(l + va) (49) 

where r is a three-dimensional Bloch vector and a the Pauli matrices. For the pur- 
pose of computation, we set the Bloch vector of the initial state to be (r sin 6, 0, r cos 9), 
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i.e. 



1 I 1 + rcosO rsini 
P(0) = i 

r sin 9 1 — r cos I 
If we subject the system to a unitary evolution 



(50) 



p(0) -► p(i) = U(t)p(0)U(tf with W(t) = exp(-|<7 3 ) 



(51) 



the Bloch vector precesses about the z-axis (at a constant polar angle of 6). For 
cyclic evolution, the parameter t takes values [0, 2tt}. Now, the matrix (46) can be 
diagonalized: 

/ 



P(0) 



with respect to its eigenvectors: 



II >= 



1 




(52) 





A 




( 




cost 






sin^ j 






|2 >= 






sin 6 


) 




{ 


— cos 9 J 



(53) 



In (52) we assume < r < 1. The geometric phase can obtained by evaluating (24): 



7g[W] = arg 



1 + r 



arctan(r tan 



e 



1 — r 



(54) 
(55) 



where f2 = 2-7r(l — cos 9) is the solid angle subtended by the Bloch vector. It 
is interesting to note that the same expression has been furnished in [10] where 
parallel transport conditions were employed. Here we have obtained the geometric 
phase without the need for such conditions. 

We now turn to the degenerate case. Here we choose a three-state density matrix 
given by: 



p= o(i+E r ^) 
6 i=i 



(56) 
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where {r{\ are the components of an eight dimensional vector and {\{\ are the 
generators of SU(3): 








1 










— i 







1 
















1 


Ai = 


1 








A 2 = 


i 








A 3 = 





-1 





A 4 = 














































1 

















— i 


























i 








A 5 = 











A 6 = 








1 


A 7 = 








—i 





1 







i 













1 










i 











-2 



For definiteness we choose the initial state to be 



P(0) 



oj 
u 
1 - 2lo 



With < UJ < 



(57) 



Here the degeneracy is manifest with the little group given by U(2) x U(l). It is 
interesting to note that in the above basis for the S77(3) generators, the subgroup 
577(2) x U (1) is spanned by {Ai, A2, A3, A 8 }. 

For an explicit evaluation of the geometric phase we consider a unitary matrix 
of 5*7(3) given by 



U(t) 



-itx 



where X = aAg + 6A4 



(58) 



where a and b are arbitrary real parameters. In a matrix form, we have 



U(t) = e^Vs 



c cos 



ct 



— iy/3i 



a sin 



ct 





iV3a t 



^sin^ 

a 2 



^sin# 

a 2 



e 2 






c cos I + i\^3a sin ^ 



(59) 
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where c = yj 3a 2 + 46 2 . The cyclicity conditions can be obtained by restricting t to 
lie within the interval [0, r] where r is evaluated by requiring that rj(r) = rj(0) with 



n(t) = -Tr[p(t)A,] 

= ^Tr[W(t)p(0)W(t)tA l 



For the U matrix considered above, we have 

2im 



(60) 
(61) 

(62) 



In the following we restrict ourselves to n = 1. With this established, the T matrix 
can be evaluated and it takes the form 



T[U- t] = e 



c cos ^ + iy/3a sin ^ 












ia\/3 t 

e 2 










-ia\f%T 



(63) 



with which the geometric phase follows by a direct evaluation of (30): 



7g[W] = arctan 



( 



sin arctan 



• tan i 



) 



2lu 



2uj 



^j- + cos(arctan[A:tan0]) 



where 



k = 



y/3a 



and 



= 



vr - 2c 



(64) 



(65) 



c cv^ 
In passing, we like to emphasize that a gauge transformation corresponding to any 
element of the little group V £ (7(2) x U(l) leaves the geometric phase invariant. 
For instance, if we had choosen 



U\t) = U{t) ■ e~ idXlt 



(66) 



which corresponds to a gauge transformed unitary matrix, then the T matrix takes 
the form 



F[U',t] 



e 2v^ 



(ay/3isin ^ + ccos ^) cos(dr) (ay/3 sin ^ — iccos ^) sin(dr) 



iaTs/3 . . . 

—ice 2 sm(ar) 



2 2 

ce 2 cos(ar) 



a 
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Evaluation of the geometric phase yields expression (64), thereby demonstrating the 
gauge invariance of (30). 

5 Concluding Remarks 

To summarize briefly, we have considered the issue of phase holonomy of both non- 
degenerate and degenerate mixed states undergoing unitary evolution. Starting with 
the non-degenerate case, we have shown that the usual procedure of subtracting the 
dynamical phase from the total phase to yield the geometric phase, does not hold for 
mixed states. To this end, we have furnished an expression for the geometric phase 
that is gauge invariant. In other words, unlike the procedure given in Ref.flO], we 
do not have to invoke the parallelity conditions to obtain the geometric phase from 
the total phase. In fact an obvious statement of gauge invariance is that any unitary 
matrix that realizes a path for the density matrix will yield a value that depends 
only on the path traversed in the state space. Moreover the form presented here also 
exhibits the parallelity conditions clearly. In essence, we have provided an explicit 
justification to the sufficient conditions furnished in Ref.[10]. The analysis is also 
extended to cover states that may exhibit degeneracies. Here with the holonomy 
taking on a non-abelian character, we have provided an expression for the geomet- 
ric phase that is manifestly gauge invariant. As in the case of the non-degenerate 
case, the form also displays the parallelity conditions which are shown to be gener- 
alizations of the ones provided in [10]. We have also furnished explicit examples of 
the geometric phases for both the non-degenerate as well as degenerate mixed states. 
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6 Appendix 



In this appendix we show that for a gauge transformation of the form (41), the 
expression for the geometric phase (30) remains unchanged. We proceed to prove 
this by establishing a few lemmas. 

Lemma la 

Let Xi[U;t] and X 2 [U;t] be submatrices of PoU(t) associated with the subspaces 
7i n and Ti m respectively: 

Xi[U;t] uu = Ip U(t)] ui/ =< n\pdA(r)\v > p,v=l,2,...n 
X 2 [U;T] ij = [poW(r)]y =< i\poU{r)\j > i, j = n + 1, n + 2, ...N 

then 

Ty h { Po U(t)F[U-t]) = Tr^X^T^J^T]) + Tr Hm {X 2 [U;T]T Hm [U;r]). 
Proof: 

The result follows immediately by taking 

n N 

Tv h { P vU{t)F[U;t}) = ]T < p\p U(r)^[U; r]\p > + ]T < j\poU(r)^[U;r]\j > 

[i=l j=n+l 

and noting that T[U;t] is block diagonal (see (39)). 
Lemma lb 

Under a gauge transformation (41), the submatrices Xi[U; r] and X 2 [U; t] transform 
as follows: 

X 1 [U;t]^X 1 [U';t] = X 1 [U;t]V 1 (t) 

X 2 [U;t]^X 2 [U';t] = X 2 [U; r]diag(e i ^+ 1 ( T ), e ie "+ 2 e id ^) 

Proof: 

This again follows from the block diagonal form of V(t) 
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Lemma 2a 

Under a gauge transformation (41), jF^ m [^/;r] transforms as 

F h JU;t] - ThJA'\t\ = diag(e^+ 1 M, e ^+ 2 M,...,e^«)t^ m [ W;r ] 

Proof: 

In proving this result, we will begin with (36). Now under U(t) — ► W(t) = U(t)V(t) 
we have 

= i[V t (^(t) t (^(t)V(t) + U(t)V{t))] kk 

= i[U{t) ] U{t)] kk -9 k k = n + l,n + 2, ...N. 

This in turn implies that 

- > 

from which we surmise that 

P'k = Pk — Ok + £,k 

where £ k is an arbitrary t-independent term. This term, however, vanishes when we 
impose initial conditions (3 k (0) = (3' k (0) = together with 0^(0) = 0. Thus, we have 
for the T- matrix 

FhJA'\t] = diag(e^+iM,e i/3 "+2M,...,e^W) 

= diag(e ifl " +1 M, e tf »+ 8 «, e^^)^ [W; r] 

which completes the proof. 

Lemma 2b 

Under a gauge transformation (41), ^ n [W;r] transforms as 

^„[W;r] - ^ n [W , ;r]=V 1 (r)t^ n [W;r] 
20 



Proof: 

We first note that under U{t) -> U'{t) = U(t)V(t), A = U(tfU(t) transforms as 

A A' = U'(tfti'(t) = VUV + V f V. 

Restricting A to the subspace H n , the n x n sub-matrix A with components 
^4 M1/ =< /i|»4|^ > transforms as 

A^ A! = v\AVi + v\Vi. 

Then from (35), with U -»• W, 

a = — „4 a 

= -{v\AV 1 + v\Vi)a' 
= -V\{AV 1 + V 1 )a' 

and thus 

± [ Vl a f ] = -Ma'. 

As a consequence of (35), we note that 

Vxq' = aV =► a'(t) = v}(t)a(t)V 

where Vo is a i-independent but otherwise, arbitrary non-singular n x n matrix. 
As in the previous lemma, this arbitrariness is removed when we invoke the initial 
conditions a'(0) = a(0) = 1, together with Vi(0) = 1. As a result Vo = 1 and we 
have 

TnJU'-A = a'(r) = Vi(r)t«(r) = V^FhMt] 
which furnishes the desired result. 



21 



The gauge invariance of the geometric phase then follows immediately from the 
above lemmas: 



7 [Z/] = arg{Tr w (p W , (T)^[W / ;r])} 

= argjTr^^lWVl^JWV]) + Tr Hrn (X 2 [U' ; t]^ h JU' ; r})} 

= arglTr^X^Tp^JZ^r]) + Tr Hm (X 2 [U;T}^ Hm [U;r})} 

= arg{Tr w (po^(T)^^;r])} 

= m- 
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